Distribution of the wavefunction inside chaotic partially open systems 
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We demonstrate both theoretically and experimentally that the distribution of the wavefunction 
inside a partially open chaotic timereversal symmetric system displays significant deviations from 
the Porter Thomas distribution. We give arguments which show that this distribution resembles 
the distribution which is expected to be found in closed chaotic systems with broken time reversal 
symmetry. 

PACS numbers: 05.45. +b 

It is now generally accepted that in the balhstic regime the scattering on boundaries of a chaotic mesoscopic 
quantum dot leads to irregular scattering phenomena which are measurable during the transport of electrons through 
the dot. Usually it is assumed that ideal leads are attached to the quantum dot and the conductance G is evaluated 
using Landauer formula relating the conductance with the corresponding S matrix 

s= (;;,',) (1) 

where r, t are the reflection and transmission matrices. In terms of S matrix the conductance reads 

G=^Tr(tt+). (2) 

Essentially the same mechanism can be applied also to the transmission of microwaves through irregularly shaped 
cavities. In this case G denotes the total transmission probability (of course in this case the factor e'^/h has to be 
omitted!). The microwave experiments have the advantage that all components of the scattering matrix are obtained 
directly from the measurement, and that scattering geometries can be easily varied in a controlled manner. Therefore 
the results to be derived below will be tested using results from microwave cavities containing randomly distributed 
scatterers. 

Investigating conductance fluctuations the standard argument assumes that the S matrix belongs to some random 
matrix ensemble (usually the circular orthogonal (COE) or unitary (CUE) ensemble). There is, however, also another 
aspect of the statistical properties of the S matrix elements which relates the S matrix elements with the internal 
wavefunction. 

The S matrix maps the incoming waves into the outgoing ones. Therefore knowing S and the structure of the 
incoming part of the wave inside the leads one can easily evaluate the value of the wavefunction at the points 
where the leads couple to the resonator. Since we assume an ideal coupling of the leads on the quantum dot the 
wavefunction is smooth at the coupling points. This means that the internal wavefunction at the coupling points 
equals the wavefunction inside the leads. Therefore from a knowledge of the statistical properties of the S matrix we 
can obtain information about the properties of the internal wavefunction. 

It is the aim of this letter to develop the above heuristic arguments and to show to what extent the transport 
through the dot changes the structure of the wavefunction. Let us start with some remarks: the important role of 
the internal wavefunction during the transport through weakly open systems (resonance transport) has been used in 
[0,^ to investigate the statistical properties of the conductance. In these papers the applied strategy was however 
just opposite to what we try to do here: the conductance distribution was obtained assuming that the internal 
wavefunction is chaotic and has a Porter-Thomas distribution. On the other hand there are also works which followed 
the same strategy and used the measured S matrix to determine the internal wavefunction, see for instance It 
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has to be stressed, however, that in all these cases theinvestigations are restricted to the resonance scattering regime 
and theinfluence of the transport on the structure of the internal wavefunction hasbeen neglected. 

To begin the investigation let us first construct a simple Hamiltonian which leads to the S matrix being usually 
used as a starting point for the description of transmission fluctuations. Starting with the description of the leads we 
assume that they support M open channels and are described by one-dimensional Hamiltonians 

Hj = -^ + A, , / = 1,..,M. (3) 

Here A; is the threshold energy of the l-th. channel and x denotes the coordinate along it. Combining these operators 
into a Hamiltonian H^x for the "external" part of the system we get 

i/e.=-l^+A, (4) 

where A is a diagonal matrix describing the threshold energies of the channels, 

A = dia5i(Ai, A2, .., Am) , (5) 

and 1 is the the M x M identity matrix. The resonator is described by a Hermitian matrix Hin of size N x N, 
where N corresponds to the number of eigenenergies taken into account, with N much larger than the number of open 
channels, N ^ M. The matrix Hin is assumed to belong to the Gaussian orthogonal / unitary Ensemble (GOE/GUE) 
for chaotic cavities (see e.g. Q for these concepts). To describe the scattering we couple the resonator and leads by 
defining the Hamiltonian Hoi the whole system as 



Ui, 



where u = {ui, . . . ,um) stands for the wavefunction inside the leads and Uin describes the wavefunction within the 
resonator. u'{0) denotes the vector of derivatives of the wave functions inside the leads taken at the points of contacts 
with the resonator (i.e. at zero of each lead coordinate). A is the N x M coupling matrix. Later on we will specify the 
coupling matrix A by assuming that the coupling has a local character (point contacts). This assumption is justified 
whenever the diameter of the junction is smaller than a typical wavelength inside the resonator. 

Let us now return to the Hamiltonian H . In the form it is given by (|^) the operator is not symmetric. To make it 
symmetric an additional boundary condition is needed: 

A^u,n = -u{0). (7) 
It is not difficult to show, that under these condition the Hamiltonian _ff is a self-adjoint operator - for the proof see 
®- . 

Before proceeding further we have to specify the structure of the internal Hamiltonian matrix Hin- On the most 
general level we assume only that this matrix belongs to the Gaussian orthogonal (unitary) ensemble. Less abstractly 
it will be however helpful to construct this matrix using the knowledge of the specific properties of the resonator 
in consideration. To specify Hin we use the Hamiltonian iJ^es of the resonator (usually a two dimensional Laplace 
operator with Dirichlet boundary conditions). Let En and /,i(r) be the eigenvalues and eigenfunctions of Hres, 

Hresfair) = Enfn{r^. (8) 

Using these solutions we define a finite dimensional internal Hamiltonian acting on the space spanned by the first 
N eigenstates of Hres, 

N 

H,n^J2^-f"fn- (9) 
n=l 

The coupling operator A maps the vector u'(0) into a certain function belonging to the N dimensional internal 
space. Let u{x) — (mi(x), M2(x), . . . , UAf(x)) denote the components of the wave function in the attached open 
channels. Applying the matrix A to the incoming vector we get 

M N 

A7.'(0)= ^^u:„(0)A„„,/„ 

771 — 1 71—1 

= ai<(f)<(0) + a2<(r)M^(0) + . . . + andUryMiO), (10) 
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where aidf = ^nifr, cti, . . . , um are the coupling constants of the individual channels and df , I = 1, . . . , M 

are functions spanned by the vectors /„; n = 1, . . . ,iV. In the experiment the channels are locally coupled to the 
resonator at points ri, . . . , tm- Here ri, . . . , vm refer to the coordinate system chosen inside the resonator and 
correspond to the zero points of the coordinates inside the leads. In order to mimic this local coupling we choose the 
functions df {r), I = 1, . . . , M in a. special way which ensures their convergence to 6 {r — fi) for N ^ oo, namely 



N 

E 

i=l 



h{n)h{r) 



(11) 



In the sense of generalized functions one has limAr^oo df [r) = 5[r — n). 

At this point few remarks are necessary: First of all it has to be stressed that now the coupling matrix A is completely 
determined by M constants ai, . . . , aM- The boundary condition takes now the form < aidf,Uin >— — m;(0), 
where <, > denotes the scalar product in the iV-dimensional Hilbert space of functions inside the resonator. In the 
limit N oo this boundary conditions are thus given by the formula 



(12) 



which relates the internal wavefunction at the coupling point to the wavefunction inside the corresponding lead. If 
we assume that the contact is ideal, i.e. that the wavefunctions match smoothly at the contact, Uin{ri) = it;(0), we 
obtain a/ = — 1 for all I. The components of the coupling matrix A are then given by 



Aim — ^flil^m)- 



(13) 



Equation (|l|) becomes a starting point for a further consideration of the statistical properties of the internal 
wavefunction. 

Since the information about the structure of the internal wavefunction is hidden in the scattering data we have first 
to evaluate the corresponding S matrix. Following we solve the equation HU = EU for a scattering energy E. 
The eigenfunction 



U{E) 



u{E^ x) 



(14) 



solves the equations 



~^u{E,x)+Ku{E,x) \ ^ 
H,nuUE)+Au'{E, 0) 



E 



u{E,x) 



(15) 



For energies E > Xj, j ^ 1, . . . , M, the scattering solution inside the leads can be presented in the form 



uiE, x) 



„i\/E—Ax 



</E^ 



</E^ 



(16) 



where AincAout are the amplitudes of the incoming and outgoing wave, respectively. For all E > maxJXj} every 
solution is bounded and the scattering matrix S{E) can easily be defined. The normalization used in ( |l6[) ensures 
that the S matrix relates the amplitudes of the incoming and outgoing waves as 



The scattering matrix can be calculated substituting ( [l^ ) into ( |l5| ) and applying the boundary condition 
leads to 



SiE) 



W^{E-H„ 



W^E-H^, 



(17) 
which 

(18) 



where W — AQ ^ and Q = Q{E) denotes the M x M matrix Q{E) = {E — A) 4 . This S matrix can be rewritten as 



S{E) = 1 - 2iW'i 



E ~ H, 



-W, 



eff 



(19) 
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where -ffe// = Hm+iWW^ stands for the effective, nonhermitian Hamiltonian. The assumption on the point character 



of the contact leads in the Hmit ^ cxd to Heff = Hres + *Sn=i Q^nV^ ~ 5{r — tVi). Effective Hamiltonians of 
this type has been used previously on a heuristic level in for the description of conductance fluctuations in quantum 
dots. In the present case, however, the coupling of the imaginary effective potential depends on energy (frequency) 
of the incident wave. This fact has been observed experimentally in 0|. 

For the value of the external component u{E, x) at the coupling points we get from (|l( 



u(£;,o) = Q(i-5M™c; (20) 

and finally from (|l2|) a relation between the S matrix and the internal function: 



Uin{r2) 



\ Um(nVf) / 



Q{S-l)Anc. (21) 



The relation ( pi| ) is valid not only for the case of a resonant scattering but also in situations with many overlapping 
resonances. This can be understood when the formula is compared with previous results, where a relation between 
the resonance wavefunction and the S matrix was obtained for the case of one isolated resonance • The relation 

( pT| ) gives us the possibility to look for the structure of the internal wavefunction in situations where the S matrix 
belongs to certain ensembles of random unitary matrices. 

In a typical experiment the wave is fed into the resonator through one of the attached channels - say through the 
first one - and the wavefunction is measured either by measuring the reflection in the entrance lead or the transmission 
to an exit lead. This means that Ainc — (1,0, ..,0) and 

u^n{rk) = j^-;^^j-y^ {Sik - Sik) (22) 

with fc = 1 in the reflection case and fc > 1 in the transmission case. Assuming At to be constant (i.e. not varying 
from sample to sample) we obtain finally that the statistics of Ui„(rfe) coincides with the statistics of the S matrix 
elements Sik- 

It is widely accepted that the structure of the S matrix of a chaotic system with time reversal symmetry is described 
by a generalized circular ensemble (so called Poisson kernel) introduced by Gaudin and Mello jsj, see also @. The 
structure of this ensemble is fully determined by the mean value < S* > of the S matrix. The mean value describes 
the "quality" of the contact. It is related to the coupling matrix W by [0 

< Su >= (23) 
1 + 7i 

with 7; given by 

N 

(24) 

n=l 

Note that the orthogonality relation (|^ follows from (|3|) and is satisfied whenever the distance between the 
coupling points exceeds the typical wavelength in the scattering problem. In the case of ideal coupling we have 
< S >— and the ensemble coincides with the Circular Orthogonal Ensemble of Dyson. The relation between the 
Poisson kernel and S matrices of the form (|l^) is described in 111. 



The statistics of the corresponding S matrix elements Ski was evaluated in |12|. It was shown that for sufhciently 
large M the S-matrix elements can be regarded as statistically independent. Moreover the real and imaginary part 
of Ski are independent Gaussian distributed variables. 

Let us first consider a measurement of the transmission between leads 1 and 2. Denoting 5i2 = X + iY and following 
[O we obtain 



with a, b given by 



P^S,,) = ^e-^^'e-"^' (25) 
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1 + <S11><S22> 

^ (26) 



~ (l-<Sll>^)(l-<S22>^) 

Ij _ 1-<511><522> 



(l-<Sll>2)(l~<S22>=^) 

and < 5*11 >, < S'22 > being the mean values of the corresponding matrix elements (they are assumed to be real for 
simplicity) . Finally using (E3) we obtain for the distribution of the modulus of the normalized internal wavefunction 
at point r2 IH] 

P (|m„P) = Zexp (-Z2|u,„|2) /g (^ZVZ^ -l\u,n\^) (27) 

with 

Z = , 28) 

VI- < ^11 >'< S22 

A similar relation holds also for the reflection measurement. The internal wavefunction is studied using the matrix 
element 5*11 which measures the direct reflection of the wave into the incoming channel. For the distribution of 
Sii ^X + iY one has pi: 



with the coefficients a, b given by 



P^Sn) = ^e-''(^-<^»'e-^^' (29) 



„ _ 1+<S11>^ 

"-Ji-<sii>2)^ (30) 

and with <S'ii>=<X> denoting the mean of the S matrix element. Using that we get for the distribution of the 
normalized internal wavefunction 

< M„ >= 1- < Sii > (31) 
and for the distribution of Ui„ = Uin~ < Ui„ > 

P (|m„|2) = Zexp {-Z^\u„f) lo [ZVZ^ - l|u„,p) (32) 



with Z equal to 



Z = =. (33) 

VI- < ^11 



In the COE case with < 5 >= (ideal coupling of channels) both distributions (^ and ( p2| ) lead to Poissonian 
distribution: 

P((|uz„|') =exp(-|u,„|2). (34) 

In the other extreme case of very weakly coupled channels for which we have < >~< S22 >~ ±1 and hence 
Z 00, which reduces ( p7| ) to the Porter-Thomas-distribution 

P(|.„P) = _l=exp f-^^) . (35) 



2 

in I 



The same holds also for the distribution ( p^ ) with < Sw >~ ±1. 

It is worth stressing that in the COE case the distribution of the internal wavefunction coincides with the distribution 
of eigenvectors of an GUE ensemble, i.e. with a case of fully broken time reversal symmetry. 

For smaller number of channels the distribution of the S'-matrix elements does not factorize and the distribution of 
the internal wavefunction cannot be found explicitly. Nevertheless it becomes clear that it depends strongly on the 
number of open channels. For instance in the case of two ideally coupled channels (COE case) we obtain from |Q 

^(i"™n = T7^^- (36) 
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Microwave billiards are especially well suited to test the predictions on distributions of the S'-matrix elements 
(seeEqs. (^ ) and (p9|)) and of transmission and reflectionprobabilities (see Eqs. (|27| ) and (^2|)). AUquantities entering 
into these expressions are directly available from theexperiment. It is of special importance that there are only two 
parameters, namely the averages < Sn > and < S22 >■ As these quantities, too, canbe measured, there is no free 
adjustable parameter in the theory. 

The measurements were performed in a rectangular microwave resonator with 10 to 17 randomly distributed scat- 
terers; fig. |^ shows a typical arrangement. The number of scatterers should be sufficient to block most of the bouncing 
ball modes in order to make the system chaotic; on the other hand the mean distance should be at least of the order 
of the typical wavelength. The measuring technique was described earlier 
360B, Wiltron company, was used supplying real and imaginary part of all components of the S matrix. 

To study the transport through the resonator 16 antennae (thin copper wires of diameter 0.2 mm) were put into the 
resonator. The antennae act as single scattering channels as their diameter is small compared to the wavelength in the 
total frequency range. Only two antennae were really used for reflection and transmission measurements. Through 
the antenna 1 microwaves are coupled to the resonator and reflection (511) is measured. The transmission (512) is 
measured with the help of the antenna 2 . All other antennae not used are closed by 50 H. loads and act as drains for 
the microwaves. 

In Fig. 1^ parts of a typical measured S'n-spectrum are shown. The qualitative difference between the non- 
overlapping regime at lower frequencies (a) and the overlapping regime at higher frequencies (b) is immediately 
evident. For the case of non-overlapping resonances we performed about 600 measurements with different positions 
of ten scatterers from 2 to 4 GHz with a resolution of 1 MHz. In the region of well separated resonances one can 
derive from the billiard equivalent of the Breit-Wigner formula a relation between |ui„p and the measured 5*1 1 at the 
maximum of the resonances, 

|u,„|2cxr(l-Re(S'ii)), (37) 

where F is the width of the resonance jj]. We fitted the single resonances in the real part of Sn with a Lorentzian 
and extracted height and width of them to obtain jiimp. In total we obtained a sample of more than 12000 values 
for I Mm p. The histogram of these values is shown in Fig. ^ and compared to the theoretical curves for GOE and 
GUE and the distribution ( |3^ ) using Z « 1.17, calculated from the experimentally obtained < >~ 0.72. The 

experimental distribution fits very well to the theoretical one. It should be noted that no parameter had to be fitted 
to obtain this accordance. 

At higher frequencies where single resonances are no longer well resolved we can take the total spectrum for the 
determination of the S-parameters. The calibration poses a problem. By application of standard procedures the 
influence of cables, connectors etc. is efficiently calibrated away. The influence of the antenna wire itself, however, 
cannot be removed by the calibration and results in a long range variation of baseline and phase over several GHz. 
Apart from some regimes, which were excluded from the further analysis, the drift of the phase could be corrected away 
by a polynomial background subtraction. It was not possible, however, to discriminate between phase drifts from the 
antenna and real phase shifts by the billiard. Therefore the determination of the phase of 5*11 is necessarily erroneous. 
This can be seen in Fig. ^,b where the distributions of Re(S'ii— < >) (a) and of Im(5ii) (b) are plotted. The 
error in the phase determination is responsible for the clear deviation of the distribution of Re(5'ii— < Sn >) from 
a Gaussian. For Im(S'ii), on the other hand, exactly the expected Gaussian behavior is found though here of course, 
too, the phase determination is incorrect. But as the average of Im(S'ii) vanishes in contrast to that of Re(S'ii), an 
error in the phase determination is not able to disturb the Gaussian behaviour. It should be noted that errors in the 
phase do influence only the distribution of Re(5ii— < 5*11 >), for all other distributions discussed here this phase is 
not of relevance. 

In Fig. ^ the histograms obtained from the measured spectra for the distribution of the internal wavcfunction |Mi„p 
for (a) Sii (reflection measurement) and (b) 5*12 (transmission measurement) in the region of overlapping resonances 
are shown. The results are compared with the theoretical predictions ( ^2|) and ( ^7| ) with the parameter Z « 1.06 
calculated from the mean value < Sn >. As we did not measure ^22 we assumed < Sn >~< S22 > to calculate Z 
from Eq. (^). Since the geometries of the antennae are identical, this assumption seems to be plausible. In Fig. Uc,d 
we show the distributions of Re(5'i2) (c) and Im(S'i2) (d). Again the Gaussian shapes are found in accordance with 
the theory. 

The experiments have shown that the ansatz of random unitary matrices by Pereyra and Mello [ p^ for the scattering 
matrix can perfectly account for the observed distribution of S'-matrix elements found in our microwave billiard 
(apart from one point where imperfections in the calibration make the comparison impossible). The change in the Z 
parameter from Z = 1.17 in the case of separated resonances to Z — 1.06 in the case of overlapping ones shows further 
that the presence of the 14 not used antennae closed by 50 fl transforms the system from the GOE to, essentially, 
the GUE behaviour. In the region of overlapping resonances we found transmission and reflection behaviour as being 



15 3]. A vector network analyzer, model 
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already essentially a GUE-like (see Fig. ||). In other words: the presence of the drains in the form of closed antennae 
transforms the standing waves of the original billiard, more or less completely, into running waves propagating from 
the entrance antenna to the different exit ports. 
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FIG. 1. Rectangular billiard (45 cm * 20 cm, height 0.8 cm) with 13 to 17 movable scatterers of diameter 2 cm positioned at 
random on a 18*8-grid. There are 16 fixed antennas (copper wires of diameter 0.2 mm) connected to the billiard at randomly 
chosen points with a minimum distance of 2.5 cm. 



FIG. 2. Parts of a typical S\i spectrum after calibrating away the effects of cables and connectors. At lower frequencies the 
resonances are sharply separated whereas at higher frequencies they overlap. 



FIG. 3. Histogram of the |ui„ |^-values at low frequencies (non overlapping regime) in a semilogarithmic plot. The dashed and 
dotted lines correspond to Poisson distribution [Z — 1, GUE) and Porter- Thomas distribution {Z — > oo, GOE), respectively. 
The solid line depicts the distribution (p2|) with the parameter Z ~ 1.17 calculated from < >. 



FIG. 4. The distributions of real and imaginary parts of Su — < Su > (a,b) and S12 (c,d). All distributions are normalized 
to have the variance equal to one. The solid lines correspond to the expected Gaussian behaviour. 



FIG. 5. Histograms of |winp for (a) reflection (Su) and (b) transmission (512) in a semilogarithmic plot. The dashed and 
dotted lines correspond to Poisson distribution {Z = 1, GUE) and Porter- Thomas distribution {Z 00, GOE), respectively. 
The solid line gives the distribution (p3) with the parameter Z ~ 1.06 calculated from < 5ii >. 
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